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ABSTRACT 

Numerical  calculations  are  carried  out  on  a large  number  of  aromatic-  and  polyene- 
C/"%ype  conjugated  molecules  by  the  methods  of  the  flrat  paper  of  this  series.  These  cal- 
' ^“^culations  include  all  energy  levels,  the  wavelengths  of  the  first  four  electronic  tran- 
sitions, the  complete  w-electron  wave functions,  the  bond  densities  and  bond  orders, 
free  valence  and  inactivities,  and  dipole  moments.  A method  of  estimation  of  bond 
lengths  with  a root  mean  sqxiare  error  of  O.OI6X  is  presented.  Wherever  possible  com- 
parisons are  made  with  experiment  and  with  the  results  of  the  LCAO  molecular  orbital 
theory  in  the  approximation  developed  by  Coulson  and  Longuet -Higgins . The  results  in 
general  are  at  least  as  good  as  those  of  the  LCAO  calculations. 


THE  FEIST  PAPER  of  the  present  series  has  presented  a detailed  consistent  mathematical 


frameworic  for  the  application  of  the  free-electron  model  to  conjugated  ir-electron 
systems.  The  present  paper  presents  the  results  of  a series  of  calculations  designed 
to  demonstrate  the  applicability  and  test  the  potentialities  of  the  free-electron 
model. 


ftAny  of  the  previous  theoretical  calculations  on  conjugated  hydrocarbons  have  been 
done  by  the  LCAO  (linear  combination  of  atomic  orbitals)  molecular  orbital  theory  in 

P i 

the  approximation  developed  by  Coulson  and  Longuet-Hlggins . It  is  emphasized  that 
this  method  will  be  referred  to  throughout  this  paper  as  the  LCAO  method,  and  that  this 
abbreviation  will  refer  only  to  the  Coulson-Longuet-Hlgglns  approximation  Just  men- 
tioned. It  has  been  pointed  out  in  I that  the  LCAO  three-dimensional,  one-electron 
model  actually  closely  approximates  a one-dimensional  model,  such  as  the  present 
theory,  so  it  is  not  surprising  that  there  exists  a very  close  equivalence  between 

LCAO  results  and  the  results  of  the  present  paper.  For  this  reason  a detailed 
^Thls  work  was  assisted  In  part  by  the  Office  of  Naval  Research  under  Task  Order  IX  of 
Contract  K6ori-20  with  The  University  of  Chicago. 

^K.  Ruedenberg  and  C.  W.  Scherr,  J.  Chem.  Phys.  21,  000  (1953),  hereinafter  referred 

to  as  I. 

^C.  A.  Coulson  and  H.  C.  Longuet-Hlggins,  Proc . Roy.  Soc.  (London)  39  (1947). 
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comparison  of  the  two  theories  is  presented  at  the  end  of  this  report. 

The  complete  rr-ele<!tron  wavefunctions  for  each  of  the  molecules  considered  here 
are  presented.  These  wavefunctions  will  be  useful  to  calculate  the  7r-electron  density 
at  any  point i to  make  perturbation  treatments,  and  to  calculate  oscillator  strengths  of 
the  electronic  transitions. 


In  referring  to  particular  atoms,  the  standard  organic  numbering,  as  indicated 
on  Pig,  1,  part  (a),  is  followed}  bonds  are  referred  to  by  their  terminal  atoms.  For 
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Pig.  1.  Bond  densities  and  bound  valences.  The  members  of  each  pair  marked  (a)  pre- 
sent the  bound  valences  of  the  atoms  indicated.  The  members  marked  (b)  present  the 
total  mobile  bond  densities.  The  bound  valences  for  joint  atoms  of  the  altei’nant  mole- 
cules are  not  presented}  they  are  generally  greater  than  2.0,  usually  around  2.J. 


example  the  bond  connecting  the  a and  p positions  in  naphthalene  is  referred  to  as 
bond  1—2. 
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1.  THE  SECULAR  EQUATION  AND  THE  ENERGY  LEVELS 
No  particular  difficulty  was  encoimtered  In  solving  the  secular  equations  for  the 
eigenvalues . The  Ijirgest  determinant  after  factorization  according  to  symmetry  species 
was  9x9  (ovalene).  The  first  n energy  levels  (n  Is  the  number  of  Tr-electrons)  and 
their  symmetries  are  presented  In  Table  I.  The  remaining  levels  and  their  symmetries 
are  trivially  found  from  these  tabulated  values.  By  the  pairing  theorem  for  energy 
levels  of  alternant  molecules  (see  paper  1),  for  every  solution  where 

Kj  - 2ffD/X^  , 

D «■  the  internuolear  distance, 
and 

Xj  ■ the  wavelength  of  <i^,  the  1th  7r-eleotron  wavefuhotion,  . 

there  exist  solutions  where  tn  may  take  any  Integral  value,  positive  or 

negative,  or  zero  and  with  the  restriction  that  > 0>  The  d^th  level  between  mn 

and  (m+l)ff  (altogether  the  (mn+l)th  level)  has  the  same  symmetry  as  the  1th  level 
between  zero  and  rrj  the  energy  relation  being,  as  indicated  above,  ■ mTT+fCj^. 

The  energy  levels  of  nonelternant  molecules  are  not  mirrored  at  half-integral  multiples 
of  rr)  that  is,  for  every  there  only  exist  solutions  ••  mn+Xj^,  m ■ 0,  *1, 

Instead  of  k^,  the  value  aotvutlly  entered  In  Table  I is  oosk^.  The  symmetries 
tabulated  refer  to  the  axes  of  Fig.  2 and  az<e  convenient  for  purposes  of  calculation. 


Pig.  2.  The  arbitrary  coordinates. 
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■3 

The  group  theory  notation  la  also  referred  to  these  same  axes.'^ 

2 . THE  WAVEFUNCTIONS 

The  wavefunotlon  of  the  _lth  orbital,  may  be  represented  by 

dj  n £ , (2) 

B 

where  la  that  part  of  located  on  branch  B and  la  zero  elsewhere.  Each 
then,  is  of  the  form 

*B1  " aBi°08()Cj^x+6gj^) , (2a) 

where  is  the  phase  of  i<)gj  referred  to  an  (arbitrary)  coordinate  axis  located  on 
branch  B. 

The  complete  m-eleotron  wavef unctions  of  all  moleoulea  considered  In  this  paper 
are  presented  In  Tables  1,  II,  and  III.  Table  I lists  oosk^.  Table  II  lists  agj^®,  and 
Table  III  lists  tanfig^.  Figure  2 depicts  the  numbering  of  the  branches  and  the 
origins  of  the  Xg's.  The  bases  of'the  arrows  are  the  origins,  and  the  heads  indicate 
the  direction  of  positive  Xg.  In  no  calculations,  except  those  of  the  spectral  lines, 
was  It  necessary  to  know  or  6g^  explicitly.  Wherever  the  phase  is  obvious  by  a 
symmetry  consideration,  no  entry  has  been  made  In  Table  III. 

There  are  two  methods,  described  in  the  Appendix,  which  may  be  employed  to  find 

p 

the  wavefimotions . The  first  method  consists  of  solving  for  Sg^  and  tanfig^  directly 
from  the  Joint  conditions  [I,  Eqs.  (1.32),  (1.33),  and  (1.3^^)]  and  from  Eq.  (2.7)  of  I. 
The  second  method  consists  of  finding  the  eigenvectors  of  the  secular  equations  by  the 
usual  methods,  [see  I,  particularly  Eqs.  (2.36)  and  (2.37)]  as  In  LCAO  calculations. 

The  first  method  Is  generally  most  convenient  If  the  molecule  possesses  symmetry.  The 
second  is  generally  more  convenient  If  the  molecule  Is  large  and/or  possesses  little 
or  no  symmetry  The  only  molecule  for  which  the  second  method  was  used  In  the  present 
calculations  Is  ovalene. 

The  second  method  only  gives  the  value  of  the  wavefunotlon  at  the  atoms.  The 
n-Blectron  density  at  Che  mid-point  of  a bond  is  found  from: 

^The  Z axis  Is  always  taken  perpendicular  to  and  out  of  the  plane  of  the  paper  except 
for  molecules  with  Cg^  symmetry,  where  the  Z axis  is  as  indicated  in  Fig.  2 and  the 
X axis  Is  perpendicular  to  the  plane  of  the  paper.  The  orbitals  are  antisymmetric  with 
respect  to  6^,  designating  reflection  at  the  plane  of  the  molecule. 

**The  phase  angle  is  arbitrary  by  an  additive  factor  imm,  m Integral.  One  takes 
ir/2  ^ 6g^  > -tt/2. 
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value  Ox  (p.  at  bond  «•  . 

^2  COS^Kj 

The  quantities  of  Tables  II  and  III  are  then  found  from 

Bi  «/_4__  2/ 


tani. 


2(ain)Cj)iJigj  (0) 


and 


« (l+tan®6gj^)(>gj_^(0)  , 


(5a) 


(5b) 


(3o) 


where  0j^(r)  is  the  value  of  the  ith  orbital  at  atom  r,  for  example,  and  in  (3b)  and 
(3o)  Pj^{x)  la  itfl  value  in  branch  B at  coordinate  x. 

3.  SPECTRA 

The  wavelengths  of  radiation  for  the  four  tranaitions  between  the  two  highest 
filled  and  the  two  lowest  unfilled  energy  levels  are  presented  in  Table  IV.  The 
nomenclature  and  assignment  is  due  to  Platt  and  Klevens,^  These  wavelengths  are  oom> 
pared  with  the  center  of  gravity  of  the  corresponding  Singlets  and  triplets.^  Unfor- 
tunately, most  of  these  triplet  data  are  estimates.  Where  the  estimates  of  the  trip- 
lets were  not  available,  only  the  singlets  are  listed.  There  will,  of  course,  always 
be  triplets,  which  generally  lie  in  such  a direction  as  to  Improve  the  agreement  of  the 
theoretical  calculations.  Some  corresponding  LCAO  calculations^  are  also  tabulated  for 

comparison. 

The  free-eleotron  model  (and  the  LGAO  model)  gives  agreement  to  within 
10,000om~^,  usually  better.  Indeed  the  prediction  of  the  longest  wavelength  transi- 
tion is  generally  quite  good.  These  results  are  the  more  remarkable  when  one  considers 
that  there  is  only  one  (o.uasi-)  adjustable  parameter  in  the  free-eleotron  model,  D, 
the  internuclear  distance.  Actually,  no  attempt  at  adjustment  was  made]  D has  been 
taken  to  be  1.39^  (the  benzene  C-C  distance)  in  the  aromatic  systems,  and,  1.34^  (the 
ethylene  C-C  distance)  in  the  n-polyenes.  When  comparing  the  LCAO  model  predictions 
with  the  free-eleotron  model  predictions,  it  should  be  borne  in  mind  that  in  the 
former  there  are  two  adjustable  parameters,  S and  P;  S is  usually  fixed  at  about  0.25, 
and  p is  varied  to  give  a best  fit. 


B.  Klevena  and  J.  R.  Platt,  J.  Chem.  Phys.  1^,  470  (19*t9)j  J.  R.  Platt,  J.  Chem. 
Phya . 18,  1168  (1950). 

^K.  Ruedenberg  and  J.  R.  Platt,  private  communication. 
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4,  BOND  DENSITIES 

The  total  mobile  bond  order,  defined  by  Coulson  and  Longuet-Hlgglna^  for  their 
LCAO  theory  would  be,  In  free-eleotron  terms,  as  follows! 

n/2 


Ts 


2 £ d.(r)0,(s) 
J-1  ^ ^ 


(4) 


where  Is  the  total  mobile  bond  order  of  bond  r-e,  refers  to  the  orbital  and  0j(r) 
and  0j(s)  are  the  values  of  the  J[th  orbital  wavefimotlon  at  atoms  r and  s_,  respectively. 
This  definition  Is  an  approximation  to  the  electron  density  In  the  bond. 

The  total  mobile  bond  density  for  the  free-eleotron  model  is  defined  by 

n/2 


b(rs)  ■>  2 £ (r-e) 


J-1 


(5) 


(6) 


where  ^j(r-e)  Is  the  value  of  the  .J_th  orbital  at  the  midpoint  of  bond  r-s.  This  defi- 
nition is  also  an  approximation  to  the  electron  density  in  the  bond  which,  however,  is 
consistent  with  the  definition  of  the  electron  density  at  the  atoms.  For  comparison 
with  Eq.  (4),  Eq.  (5)  may  be  written  approximately  as 

n/2 

b(rs)  » f[ft  £ 0,(r)^,(s)]  , 

J-1  •’  ^ 

where  f is  some  constant,  equivalent  to  a mean  value  of  1/oos^^r^,  and  may  be  taken  as 
« 1.25.  It  would  seem  that  the  most  exact,  but  probably  unnecessary,  definition  for 
the  electron  density  would  be  an  Integrated  density  extending  half  a bond  length  to 
either  side  of  the  point  In  question. 

The  bond  densities  calculated  from  Eq.  (5)  are  presented  In  Figs.  1 and  ?,  part. 
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5.  BOUND  VALENCE  AND  REACTIVITIES 

CoulBon  and  Longuet-Hlgglns"^  have  defined  a useful  quantity,  the  free  valence, 
which  has  been  successfully  correlated  by  many  authors  with  chemical  reactivity.  It  is 
a second-order  effect  which  becomes  of  importance  in  alternant  hydrocarbons  because  the 
"electron  density"  at  the  atoms  (the  first-order  effect)  is  the  same  (but  v.JL. ).  The 
free  valenoe  is  the  potential  bond-forming  power  a particular  atom  in  a molecule  has. 
The  free  valenoe  is  obtained  by  forming  the  sum  of  the  bond  orders  of  the  bonds  termi- 
nating on  the  atom  under  consideration  and  then  subtracting  this  sum  from  a certain 
number,  is  obtained  in  the  following  manner:  If  one  were  to  have  available 

calculations  on  all  possible  hydrocarbons,  then  he  would  be  able  to  plok  out  the  max- 
imum possible  value  for  the  sum  just  mentioned  above.  This  maximum  value  is  taken  to 
be  Analytically  the 

free  valence  of  atom  r » N_._-  S p„„ 

max  *^rs 

all  s neighboring 
to  r 

This  concept  may  be  taken  over,  ^ toto.  into  the  free-eleotron  picture , The  present 
author  has  chosen  to  tabulate  merely  the  sums  of  the  bond  densities,  and  in  order  to 
Indloate  the  complementary 'nature  of  these  sums  to  the  free  valenoe,  has  called  them, 
following  a suggestion  of  Ruedenberg's,  the  bound  valenoe.  Hence,  the  smaller  the 
bound  valence,  the  more  reactive  the  atom  (to  eleotrophlllic  reagents).  These  bound 
valences  are  presented  in  Pigs.  1 and  3.  The  base  of  the  arrows  indicates  the  oorre- 
spondlng  atom. 

In  the  first  paragraph  of  this  section  it  has  been  mentioned  that  the  first-order 
effect  for  reactivities  is  the  "electron  density"}  however,  it  should  be  noted  that 
the  LCAO  model  has  tacitly  implied  that  the  first-order  effect  is  due  to  the  electron 
population  (see  paper  I,  since  to  the  author's  knowledge  this  distinction  between 
density' and  population  has  never  been  made  before).  This  point  is  being  investigated 
further.  In  this  paper  once  D has  been  taken  as  the  unit  of  length  there  is  no  numer- 
ical difference  between  density  and  population  for  bond  or  for  non-Jolnt  atom  den- 

■5 

sitles.  At  the  Joints  the  population  Is  ^ of  the  density.  No  bound  valences  are 
reported  for  Joint  atoms  (they  are  generally  greater  than  2.0,  usually  around  2. 3). 

In  nonalternant  molecules,  it  is  the  total  electron  density  at  the  atoms  which  Is 
reported  rather  than  the  bound  valence. 

'^C,  A.  Coulson  and  H.  C.  Longuet-Hlggins , Rev.  Sclentlflgue  85.  929  (194?). 
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Table  V compares  the  (theoretically)  most  reactive  atoms  In  each  molecule  with 
the  known  chemistry  of  the  hydrocarbons.  The  only  serious  discrepancy  with  exper- 
iment Is  found  in  the  phenanthrene  molecule,  where  the  1,8  positions  are  Indicated  to 
be  slightly  more  reactive  than  the  9,10  positions.  If  one  attempts  to  explain  this 
by  notloing  that  there  la  an  over-all  very  high  oonoentratlon  of  ir-eleotron  density 
about  the  9,10  positions  (^.e. , atoms  9 and  10  contiguous  and  connected  by  the  very 
high  density  9-10  bond)  then  a difficulty  arises  In  explaining  why  the  somewhat 
similar  pyrene  molecule  is  not  reactive  at  the  pyrene  1,2  positions.  An  exactly 
similar  difficulty  is  encountered  In  the  LCAO  model. 

6.  DIPOLE  MOMENTS 

Of  all  the  molecules  considered  here,  only  fulvene  and  asulene  can  (theoretically) 
have  dipole  momenta,  and  Indeed,  this  la  an  experimental  fact,  that  alternant  hydro- 
carbons have  at  most,  small  dipole  momenta.  The  calculated  moments  are:  Pulvene  4.4 

D (Debyes);  aaulene  6.3  D.  These  may  be  compared  with  the  (unrefined)  LCAO  values  of 

o 

4.8  D and  6.9  D,  respectively.  Experimentally,  fulvene  would  be  expected  to  have  a 
dipole  moment  of  about  1.2  D on  the  basis  of  the  dipole  moments  of  its  derivatives, 

o 

while  that  of  azulene  la  of  the  order  of  unity.  There  are  many  well-known  reasons 
for  these  discrepancies  among  which  the  C-eleotron  contribution  to  the  dipole  moment 
and  the  effect  of  the  true  bond  length  should  be  mentioned. 

7.  BOND  LEN0TH3 

CoulBon^  has  successfully  demonstrated  that  on  the  basis  of  the  total  ir-electron 
density  (Coulson's  total  mobile  bond  order)  at  the  bonds,  predictions  of  the  bond 
lengths  In  these  molecules  can  be  made  which  agree  fairly  well  with  experiment; 
usually  within  the  experimental  error.  The  present  author  has  found  the  following, 
completely  empirical  equation  to  be  satisfactory  for  use  with  the  bond  densities  cal- 
culated on  the  basis  of  the  free-eleotron  model.  It  reproduces  the  known  bond  lengths 
with  a root  mean  square  deviation  of  0.016^; 

length  of  bond  r-e(A)  = 1.665-0.1398[l+b(rs)l  . (7) 

Equation  (7)  has  been  obtained  from  the  curve  of  Pig.  4.  Figure  4 is  analogous  to  the 
well-known  figure  due  to  Coulson^  constructed  on  the  basis  of  the  LCAO  model  results, 

®LCA0  data  taken  from  B.  Pullman  and  Pullman,  Les  theories  llectronlgues  de  la 
ehlmlc  organloue  (Masson  and  Company t Paris,  1952),  pp.  399  and  406. 

^C.  A.  Coulson,  Prop . Roy.  Soc ■ (London)  A207,  91  (1951),  and  also  In  J.  Phys . Chem. 

56,  311  (1952). 
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Pig.  4.  This  curve  shows  the  relationship  between  the 
calculatei^.  total  mobile  bond  density  and  the  experimental 
data  now  available,  A,  anthracene;  Be,  benzene;  Bu, 
butadiene;  C,  ooronenej  Q,  graphite;  N,  naphthalene,  0, 
ovalenej  P,  pyrene;  Tr,  trlphenylene;  and  E,  ethylene. 


and  Indeed,  the  present  Fig. 
4 is  constructed  to  facil- 
itate comparison.  Experi- 
mental bond  lengths  are 
plotted  against  calculated 
bond  densities.  The  solid 
line  Is  a straight  line 
which  can  be  drawn  through 
the  points.^®  The  dashed 
lines  represent  O.OSoS  devi- 
ations from  the  solid  line, 
and  correspond  to  the  exper- 
imental error  for  most  of 
the  molecules  considered. 

The  present  author  has  cal- 
culated the  root  mean  square 
deviations  on  both  Pig.  4 
and  the  analogous  plot  of 
CoulBon.^^  These  deviations 


are,  for  the  free-electron  model,  O.OlSX;  LOAD,  O.OlSX.  Table  VI  presents  a brealtdown 

\p 

of  these  figures.  The  predictions  of  bond  lengths  by  the  free-electron  model  are 


^°A  better  fit  to  the  experimental  data  could  perhaps  be  obtained  if  one  disregards 
the  somewhat  doubtful  polyene  points  and  the  certainly  dubious  pyrene  points . 

Because  of  this  uncertainty  of  how  to  weight  the  experimental  data,  the  best  line 
has  been  estimated  by  an  Inspection. 

^^The  bulk  of  the  necessary  data  is  taken  from  J.  H.  Robertson,  Proo . Roy.  Soc . 

(London)  AS07 . 101  (1951).  It  has  been  necessary  to  make  a few  corrections  to  tnese 
data;  see,  for  example,  note  b to  Table  VI. 

^^The  e.perlmental  lengths  for  butadiene,  and  In  particular  the  1-2  bond  (1.46A),  are 
somewhat  dubious.  [See,  e .g. , Coulson,  and  Daudel,  and  Robertson,  Proo.  Roy . Soc. 
(London)  AgQ7,  311  (1951).]  For  a discussion  of  the  polyene  bond  lengths  from  a 
fi'cc-electron  point  of  view  see  H.  Kuhn.  J,  Chem.  Phys . 16,  840  (1948);  II98 
(1949).  See  also  Noel  S.  Bayllss,  Aust.  J.  Sol.  Rea.  M,  344  (1951).  and  earlier 
references  given  there;  his  treatment  does  not  conform  to  that  of  the  present  paper. 
For  these  and  other  reasons  the  butadiene  point  was  not  Included  In  this  cal.iUlatlon. 
With  this  point  In  they  become:  free-electron  O.OI7X;  LCAO  unchanged. 
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conslatently  better  than  the  predlotlona  by  the  LCAO  model. 

8.  DISCUSSION  OF  SPECIAL  POINTS 

Detailed  calculations  on  a series  of  n-polyenes  has  not  been  given.  One  may  con- 
veniently find  the  total  mobile  bond  density  from  the  following  formula: 

b(ra)  o 1 , (8) 

8=r+l  (n+1)  aln{TT(r+i)/(n+l)} 

where  r In  Eq.  (8)  la  the  systematic  number  of  the  rth  atom,  and  n Is  the  number  of 
atoms  In  the  polyene  chain.  This  equation  clearly  shows  the  "alternating  single  and 
double  bonds"  of  polyene  chains.  It  also  shows  that  the  strongest  bond  in  any  n- 
polyene  is  the  1-2  bond  and  Its  slater,  the  (n-l)-n  bond.  The  formula  Implies  that 
very  long  chains  would  show  only  a slight  conjugation  effect. 

In  the  oyolo-butadlene  molecule  there  exists  a filled  level,  b ^tt,  degenerate 
with  the  lowest  unfilled  level,  and,  hence,  oyolo-butadlene  has  a triplet  ground 
state. Properly  speaking,  that  Is,  there  are  two  htlf-fllled  orbitals  in  the  ground 
state,  and  Eq,  (5)  must  be  used  In  Its  general  form, 

b(rs)  = 2 n><J>/(r-e)  , (5a) 

Jnl  ^ 

iji 

where  nj  Is  the  number  of  w-electrons  In  the  J_th  orbital.  Longuet-Hlgglns  has  con- 

sidered the  following  three  hypothetical  hydrocarbons  from  the  LCAO  approach  which  also 
have,  on  that  model,  a triplet  ground  state;  l,3-(CH2)gCgH2^,  and  l,l-(CHg)gC  ■»  CHg, 
Clar's  hydrocarbon, none  of  which  have  a Kekule  structure,  but  all  of  which  have 
unstrained  0-bonds  and  an  even  number  of  Tr-eleotrons . As  might  be  expected,  all  three 
of  these  also  have  a triplet  ground  state  on  the  basis  of  the  free-eleotron  model  cal- 
culations . 

The  lowest  bound  valence  In  the  cata-condensed  ring  systems  considered  here  was 
found  In  pentacene,  which  agi’ees  with  the  observation  that  the  higher  members  of  the 
polyacene  series  are  not  only  highly  unstable,  but  progressively  more  unstable. 


^^See  for  p.'xample,  C.  A.  Coulson  and  0.  S.  Rushbrooke,  Proc.  Cambridge  Phil . Soc.  ?6, 
193  (19'*0),  p.  198. 

C.  Longuet-Hlgglns,  J.  Chem.  Phvs . l8.  265  (1950). 

^^E.  Clar,  Aromatlsche  Kohlenwasserstoffe  (Verlag.  Julius  Springer:  Berlin,  1941), 

p.  511.  This  hydrocarbon  consists  of  six  benzene  rings,  perl-condensed  Into  an  equi- 
lateral triangle. 
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9.  A COMPARISON  WITH  THE  LCAO  MODEL 

A general  statement  can  be  made  oonoerning  the  oomparleon  of  the  LCAO  model  of 
CoulBon  and  Longue t -Higgins  and  the  free-eleotron  modal  on  the  basis  of  the  oaloula- 
tions  carried  out  hero)  namely,  that  the  free-eleotron  model  itrediotlona  are  equivalent 
to  or  slightly  better  than  the  LCAO  results.  The  superiority  is  slight  and  may  not  be 


Specifically,  the  accuracy  of  the  predictions  of  spectra  seems  to  be  about  equal, 
particularly  when  the  center  of  g;ravlty  of  the  singlet  and  triplet  lines  are  considered. 
This  equality  actually  in  In  favor  of  the  free-eleotron  model,  since  (in  the  present 
calculations)  no  parameter  has  been  varied  In  order  to  give  optimal  fit,  whereas  in 
the  LCAO  calculations  at  least  one  of  the  available  parameters  has  bee-,  freely  adjusted 
to  give  optimal  fit. 

The  relative  bond  densities  In  any  series  of  related  bonds  are  both  qiialitatlvely 
and  quantitatively  reproduced  with  amazing  fidelity  in  the  two  theories,  despite  their 

O 

divergent  approaches.  Figure  5 shows  the  LCAO  bond  orders  plotted  against  the  free- 


na  tucmcN  BOODCNanr 

Pig.  5.  These  curves  indicate  the  close  correspondence  of  the 
free-electron  model  and  the  LCAO  model  calculated  bond  densities 
(orders)  for  alternant  molecules)  both  coordinates  to  the  same 
scale.  A,  anthracene;  Be,  benzene;  C,  coronene;  D,  diphenyl;  N, 
naphthalene;  0,  ovalene;  P,  pentacene;  Pe,  perylene;  Ph,  phenan- 
threne;  Py,  pyrene;  S,  stllbene*  St,  styrene;  T,  tetracene;  and 
Tr,  triphenylene. 
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electron  bond  densities,  both  to  the  same  scale.  The  points  closely  align  themselves, 
according  to  broad  generallaations  of  bond  type  (terminating  in  two  Joints,*  terminating 
in  one  joint,  etc . ) , on  lines  whose  slopes  are  very  close  to  unity.  The  deviations  of 
the  points  from  these  straight  lines  rarely  exceeds  0.02  density  units,  much  too  small 
to  be  deoisive  in  comparisons  with  experiment.  However,  there  do  exist  certain  differ- 
ences which  can  be  subjected  to  experimental  test,*  whenever  a point  on  a left-lying 
line  is  higher, than  a point  on  a right-lying  line  (of  Fig.  5),  then  the  relative  bond 
lengths  of  the  bonds  Involved  are  reversed  in  the  two  theories.  Such  pairs  of  points 
will  be  particularly  decisive  when  they  include  bonds  from  the  same  molecule.  The 
present  experimental  data  however  are  not  precise  enough  to  differentiate  and  could 
only  permit  such  comparisons  to  indicate  a trend. 

Nevertheless,  in  predictions  of  bond  lengths  baaed  on  these  electron  densities  the 
free-eleotron  model  is  superior.  In  the  important  case  of  the  naphthalene  molecule, 
important  because  of  the  high  degree  of  precision  claimed  for  the  x-ray  investigation, 
both  models  give  unsatisfactory  prediotions.  The  root  mean  square  deviations  for  this 
molecule  aret  Free-eleotron  model  0.023A;  LCAO  O.025X.  It  has  been  pointed  out  in 
Table  VI  that  the  LCAO  root  mean  square  deviation  of  O.OlsS  claimed^'  for  the 

LCAO  predictions  of  the  bond  lengths  in  this  molecule  is  obtained,  not  from  Coulson's 

0 ig 

graph,’  but  from  a special  graph  which  heavily  weights  the  naphthalene,  the  anthra- 
cene, and  the  graphite  prediotions.  These  Very  poor  root  mean  square  deviations  are 
directly  attributable  to  their  (almost  identically)  poor  prediction  of  the  length  of 
the  cross-linked  bond.^'^  As  has  been  discussed  elsewhere^*'  it  is  this  particular 
bond  which  is  subject  to  an  experimental  error  about  twice  those  of  the  other  three 
bonds  in  naphthalene. 

There  la  a possible  fruitful  field  of  comparison  of  bond  length  predictions.  As 

can  be  clearly  seen  in  Pig.  6,  the  prediotions  of  bond  lengths  in  the  odd-membered 

ring  systems  by  the  two  theories  are  widely  divergent.  Even  comparatively  rough  bond 

1 

length  measurements  on  azulene  would  suffice  for  a diacrimlnat Jon. 

. .. . ...  .. 

Coulson,  Daudel,  and  Robertson,  Prop.  Roy.  Soc.  (London)  A207.  306  (1951). 

^^The  value  predicted  by  both  models  is  about  I.IjJtX,  which  value  seems  to  be  more 
reasonable,  on  the  basis  of  a comparison  with  the  similar  bond  In  anthracene  (1.440A) 
than  the  experimental  value  (1.395^)  of  reference  11. 

18d**.  Hellbronner  (ETH,  Zurich)  has  kindly  supplied  the  author  with  some  advanced  ten- 
tative results  of  electron  diffraction  measurements  made  on  the  azuleno  molecule, 
which  indicate  that  the  cross-linked  bond  in  azulene  is  about  1.45^.  Proe-eleotron 
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OAO  OAO  OJO  080  OM  IDO  IJO  \S0  liO 
FREE  ELECTRON  BOND  DENSITY 

Pig.  6.  This  plot  shows  the  lack  of  oorresp-jndenee  of. the  free- 
eleotron  model  and  the  LCAO  model  oaioulatlons  of  the  bond  den-  . 
sitles  (orders)  for  nonalternant  molecules)  both  coordinates  to 
the  same  scale.  D,  dlfulvene)  F,  fulvens)  P,  pentalenej  Z,  asu- 
lene, 

A oomparlBon  of  the  physloal  labor  Involved  in  solving  moleoular  problems  by  the 
two  methods  is  diffioult  to  make.  The  time  Involved  for  a complete  calculation  on  one 
molecule  ran  from  one  minute  for  the  polyenes  through  four  and  a half  hours  for  penta- 
cene  to  six  days  for  ovalene. 

10.  CONCLUSION 

The  free-electron  model  has  been  shown  to  be  a versatile  and  respeotable  tool  for 
use  in  the  investigation  of  moleoular  structure.  It  is  oapable  of  extension  to  other 
systems  of  interest  than  those  considered  here  and  seems  to  be  oapable  of  extension 
to  include  electron  interaction)  it  deserves  more  attention  than  has  hitherto  been 
given  to  it.  It  seems  quite  reasonable  to  expect  that  the  theory  may  be  refined  to 
the  same  degree  as  the  LCAO  model  theory.  On  the  other  hand,  the  free-electron  model 
seems  to  have  the  same  limitations  and  inaccuracies  that  the  LCAO  method  has,  and 
further  its  present  treatment  of  polyenes  is  unsatisfactory. 
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TABLE  I 

THE  tt-ELECTRON  WAVEPUNCTIONS  t coa/c^ 

The  aymmetrlea  (SS,  SA,  eto. ) and  the  group  theory  notation  refer  to  the  axes  of 
Pig.  2.  The  aymmetrlea  Indicate  the  behavior  of  the  :irbltal  upon  reflection  at  the 
plane  of  the  molecule,  taken  In  the  order  X,  Y.  Hence,  AS  would  mean  antlaymmetrlo 
upon  reflection  at  the  XZ  plane,  aymmetrio  upon  reflection  at  the  YZ  plane;  S^A^,  aym 
metric  upon  reflection  at  any  of  three  equivalent  XZ  planea,  antisymmetric  at  any  of 
three  equivalent  YZ  planes.  A single  letter  (e.£. , in  molecules  with  symmetry) 
Indicates  the  behavior  upon  reflection  at  the  Z axis  (at  the  XZ  plane). 

Benzene  Dgyj  naphthalene  Anthracene  Vj^ 


1. 

Agu 

1.0 

1.  SS 

Blu 

2. 

AS  ' 

'"is 

0.5 

2.  SA 

®2g 

3. 

SA  j 

3.  AS 

"38 

4. 

AA 

(®2U 

-0.5 

4.  SS 

®lu 

5. 

SS  : 

5.  AA 

6. 

"28 

-1.0 

6.  SA 

"28 

7. 

"is 

-1.0 

7.  AS 

"3S 

8.  SS 

Blu 

9.  AA 

Au 

10.  AS 

"38 

Naphthacene 

Fentaoene 

1. 

SS 

Slu 

1.0 

1.  SS 

®1U 

2. 

SA 

®2g 

0.94490 

2.  SA 

"2g 

3. 

SS 

Slu 

0.79002 

3.  SS 

Blu 

4. 

AS 

"38 

0.63640 

4.  SA 

"as 

5. 

SA 

"as 

0.57008 

5.  AS 

"38 

6. 

AA 

0.53783 

6.  AA 

Au 

7. 

SS 

Slu 

0.37823 

7.  SS 

®lu 

8. 

AS 

"3S 

0.36518 

8.  AS 

"38 

9. 

AA 

Au 

0.14382 

9.  SA 

"28 

10. 

SA 

®2g 

-0.14382 

10.  AA 

Au 

11. 

SS 

Slu 

-0.36518 

11.  AS 

"38 

12. 

AS 

"3s 

-0.37823 

12.  SS 

Blu 

13. 

SA 

"28 

-0.55783 

13.  SA 

"2S 

14. 

AS 

"3s 

-0.57008 

14.  AA 

Au 

15. 

SS 

®lu 

-0.63640 

15-  SS 

Biu 

16. 

AS 

"38 

-0.7900a 

16.  AS 

"38 

17. 

AA 

Au 

-0.94490 

17.  SA 

®2g 

. 1-0 

1.  SS 

Slu 

1.0 

0.80902 

2.  SA 

"28 

0.90627 

0.59023 

3.  SS 

Slu 

0.65841 

0.42356 

4.  AS 

"38 

0.62139 

0.30902 

5.  AA 

A 

u 

0.46685 

-0.30902 

6.  SA 

"28 

0.39392 

-0.42356 

7.  AS 

"38 

0.20369 

-0.59023 

8.  SS 

®lu 

-0.20369 

-0.80902 

9.  A A 

Au 

-0.39392 

-1.0 

10.  SA 

"28 

-0.46685 

11.  SS 

Alu 

-0.62139 

12.  AS 

"sa 

-0.65841 

13.  AA 

Au 

-0.90627 

14.  a8 

-1.0 

Diphenyl 

Vh 

1.0 

1.  SS 

®lu 

1.0 

0.96385 

2.  SA 

"28 

0.92249 

0.85968 

3.  SS 

®lu 

0.57180 

0.70083 

4.  AA 

Au 

0.5 

0.64491 

5.  AS 

®3g 

0.5 

0.57601 

6.  SA 

"28 

0.31597 

0.51454 

7.  SS 

Slu 

-0.31597 

0.45582 

8.  AA 

Au 

-0.5 

0.36855 

9.  AS 

®3g 

-0.5 

0,29056 

10.  SA 

"2g 

-0.57180 

0.10681 

11.  SS 

Biu 

-0.92249 

-0.10681 

12.  SA 

"2s 

-1.0 

-0.29056 

-0.36855 

-0.45582 

-0.51454 

-0.57601 
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18.  AS 

®3S 

-1.0 

18.  SS 

®1U 

-0.64491 

19.  AA 

Ay 

-0.70083 

20.  AS 

®3S 

-0.85968 

21.  AA 

Ay 

-0.96385 

22.  AS 

®3S 

-1.0 

Stllbene 

\ 

Styrene 

^2v 

Olfulvene 

1.  SS 

Blu 

1.0 

1.  S 

®2 

0.97903 

1.  SS 

Bly 

1.0 

2.  SA 

®2S 

0.96096 

2.  S 

Ba 

0.66594 

2.  SA 

®2S 

0.89830 

3.  SS 

®lu 

0.72011 

3.  A 

As 

0.5 

3.  SS 

Bly 

0.42356 

4,  SA 

®2g 

0.5 

4.  S 

®2 

0.31308 

4.  AA 

Ay 

0.30902 

3.  AS 

®3S 

0.5 

5.  s 

=2 

-0.31308 

5.  AS 

®3S 

0.30902 

6.  AA 

Au 

0.5 

6.  A 

As 

-0.5 

6.  SA 

®2S 

0.11020 

T.  SS 

Slu 

0.24085 

7.  s 

Ba 

-0.66594 

7.  SS 

Bly 

-0.59023 

8.  SA  . 

®2g 

-0.24085 

8.  S 

Bs 

-0.97903 

8.  AS 

=3g 

-0.80902 

9.  AS 

^32 

-0.5 

9.  AA 

Ay 

-0.80902 

10.  AA 

yS 

Ay 

-0.5 

10.  SA 

®2g 

-0.84183 

11.  SS 

Biu 

-0.5 

' 

12.  SA 

®2g 

-0.72011 

13.  SS 

^lu 

-0.96096 

14.  SA 

Boa 

-1.0 

Methylene-3 

Pulvene 

°2V 

pentadlene  1,4 

®2v 

Pentalene 

^h 

1.  S 

®2 

0.96262 

1.  S 

Bg 

0.83041 

1.  SS 

Bly 

1.0 

2.  S 

0.44376 

2.  A 

Ag 

0.5 

2.  SA 

®2g 

0.70711 

3.  A 

Ag 

0.30902 

3.  s 

Bg 

0.24581 

3.  AS 

®3g 

0.43426 

4.  S 

Bg 

-0.12480 

4.  S 

Bg 

-0.24581 

4.  SS 

Bly 

0.19371 

5.  3 

^2 

-0.78158 

5.  A 

Ag 

-0.5 

5.  AA 

Ay 

0.0 

6.  A 

As 

-0.80902 

6.  S 

Bg 

-0.83041 

6.  SA 

®2g 

-0.70711 

7.  AS 

®3g 

-0.76759 

8.  SS 

®lg 

-0.86038 

Azulene 

Cgv 

Phenanthrene 

®2v 

Pyrene 

1.  S 

®2 

1.0 

1.  S 

Bs 

1.0 

1.  SS 

Blu 

1.0 

2.  S 

°2 

0.81254 

2.  A 

As 

0.89689 

2.  AS 

®3g 

0.87017 

3.  A 

0.63253 

3.  S 

Bs 

0.72015 

3.  SA 

®2g 

0.78868 

4.  S 

4 

0.39186 

4.  S 

Bg 

0.56250 

4.  SS 

Blu 

0.64306 

5.  A 

As 

0.22398 

5.  A 

Ag 

0.5 

5.  AA 

Ay 

0.5 

6.  S 

P 

-0.18536 

6,  A 

Ag 

0.33647 

6.  SS 

Blu 

0.41731 

7.  A 

As 

-0.35032 

7.  s 

Bg 

0.28404 

7.  AS 

=3S 

0.35687 

8.  S 

Bs 

-0.77636 

8.  A 

Ag 

-0.28404 

8.  SA 

=2g 

0.21132 

9.  A 

As 

-0.83952 

9.  s 

Bg 

-0.33647 

9.  AA 

Ay 

-0.21132 

10.  S 

B2 

-0.90934 

10.  s 

Bg 

-0.5 

10.  S3 

Bly 

-0.35687 

11.  A 

Ag 

-0.56250 

11.  AS 

®3g 

-0.41731 

12.  A 

Ag 

-0.72015 

12.  SA 

®2g 

-0.5 
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1.  SS 

®lu 

1.0 

2.  AS 

®3S 

0.92104 

3.  SA 

®2S 

0.82261 

4.  AA 

Au 

0.75120 

5.  SS 

Blu 

0.64460 

6.  AS 

®3g 

0.43524 

7.  SS 

Blu 

0.42469 

8.  SA 

®as 

0.42178 

9.  AS 

®3g 

0.37968 

10.  AA 

Au 

0.15986 

11.  SA 

Bgg 

-0.15986 

12.  SS 

Blu 

-0.37968 

13.  AA 

Au 

-0.42178 

14.  AS 

®3a 

-0.42469 

15.  SS 

Blu 

-0.43524 

16.  AS 

®3S 

-0.64460 

17.  SA 

Bgg 

-0.75120 

18.  AA 

Au 

-0.82261 

19  • SS 

Blu 

-0.92104 

20.  AS 

^38 

-1.0 

ovale ne 


1. 

SS 

Blu 

2. 

SA 

B2b 

3. 

AS 

=3S 

4. 

AA 

Au 

5. 

SS 

Blu 

6. 

SS 

Blu 

7. 

SA 

®2S 

8. 

AS 

B3S 

9. 

SA 

B2S 

10. 

SS 

Blu 

11. 

AS 

B5S 

1.0 

0.93890 

0.90269 

0.60139 

0.79300 

0.67428 

0.63641 

0.59895 

0.52671 

0.41907 

0.41900 
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13.  s 

B2 

-0.89689 

14.  A 

A2 

-1.0 

Trlphenylene 

°3h 

1.  S^ 

A " 

Ag 

1.0 

2.  S I 

3.  A J 

E" 

0.89610 

4.  S^ 

A2 

0.59023 

5.  S 1 

6.  A 1 

E" 

0.55400 

7.  A^ 

Al" 

0.42356 

8.  S 1 

9.  A J 

. E" 

1 

0.29074 

10.  A ] 

11.  s : 

1 E" 

-0.29074 

12.  S^ 

A * 
Ag 

-0.42356 

13 . A ' 

14.  S 

1 E" 

-0.55400 

15.  A^ 

A ” 
*1 

-0.59023 

16.  A 

17.  s 

1 E"  , 

,-0.89610 

18.  A^ 

Al" 

-1.0 

Ovalene 

12. 

AA 

Au 

0.41279 

13. 

SS 

®lu 

0.34698 

14. 

AA 

Au 

0.30085 

15. 

SA 

®2g 

0.24635 

16. 

AS 

B3g 

0.14602 

17. 

SS 

Blu 

-0.14602 

18. 

AA 

A„ 

-0.24635 

19. 

SA 

Bag 

-0.30085 

20. 

AS 

®3g 

-0.34698 

21. 

SA 

®2g 

-0.41279 

22. 

SS 

Blu 

-0.41900 

13.  AS  B.5  -0.64306 

14.  AA  Ay  -0.78868 

15.  SS  Bjy  -0.87017 

16.  AS  Bjg  -I'O 

Coronene  05^ 


1. 

S5s5 

Aau 

1.0 

2. 

AS  ) 

El- 

0.89349 

3. 

SA  J 

Ig 

4. 

5. 

SS  ■) 
AA  J 

®2u 

0.70998 

6. 

S^S^ 

Agu 

0.61693 

7. 

8. 

S^A^ 

a5s5 
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TABLE  III. 

The  reference  coorainatee  are  shown  in  Pig.  2.  Wherever  the  phase 


OrbltaU  2 

3 

4 

5 

6 

7 

8 

Molecule 

Branch 

Naphthacene 

B2 

1.0925 

0.066212 

15.631 

0.63293 

Pentacene 

B2 

-1.6333 

3.0250 

0.28906 

0.10148 

-0.25623 

-1.4330 

5.9578 

Fhenanthrene 

B1 

-0,092822 

0.41338 

-1.2236 

1.7321 

-0.99169 

1.5007 

B4 

-0.87844 

-2.8417 

1.9800 

-1.7321 

0.10223 

0.13161 

Pyrene 

B5 

-1.0895 

0.78433 

1.2107 

-0.02519 

Perylene 

B2 

0.46039 

1.5269 

-1.4366 

-0.69744 

B3 

1 .15656 

17.292 

0.81522 

0.063547 

B4 

1.3300 

-0.50211 

28.402 

1.3007 

Trlphenylene 

B2 

-6.3933 

2.1311 

-0.22518 

0.075061 

-3.6881 

B3 

-0.26968 

0.089892 

3.1611 

-1.0537 

-1.4872 

B4 

-2.5724 

0.85747 

-4.7962 

1.5987 

-0.096145 

Coronene 

B2 

0,39688 

1.2642 

B3 

1.0704 

1.2523 

B4 

-1.0704 

-1.0704 

-1.2523 

-1.2523 

2.7752 

-0.87759 

B5 

1.4055 

7.2845 

Ovalene 

B2 

0.75161 

-0.45800 

0.054668 

9.7656 

-0.38012 

-0.69565 

1.0065 

B4 

0.69905 

1.0701 

1.1320 

1.9561 

2.0230 

4.2473 

1.3591 

B5 

0.16904 

-1.1287 

-0.66845 

-0.27687  138.63 

-0.52226 

-0.39123 

B6 

-1.7143 

1.5406 

0.27330 

5.2209 

-0.039311 

1.2565 

4.7215 

B7 

5.6332 

1.3231 

3.0116 

-0.54008 

-0.73290 

2.1215 

0.33171 
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is  obvloxia  from  symmetry  oonsideratlons , no  entry  has  been  made. 


9 10  11  12 


9.7850  0.25472  -0.2350 


-0.70262 

-1.2338 

0.46804 

1.2294. 

0.49572 

0.032048 


-0.067714 

-1.4004 

-1.1245 

0.32469 

1„1245 

1.1245 

-0.32469 

-0.32469 

0.49741 

3.7931 

-0.021495 

0.38941 

-0.57068 

92.121 

-1 .8669 

0.82172 

0.29987 

1.1876 

2.1416 

-1.3251 

2.6748 

-0.27299 

-0.44961 

-1.7400 

0.35086 

2.2731 

1.7960 

-1.2067 

0.09435 

1.5511 

13  14  15  16 


0.16745 

-0.21940 

1.4730 

■0.35993 

1.2520 

-0.29737 

0.049343 

0.64402 

-1.7090 

1.0954 

-0.25094 

3.2016 

-1.5176 

0.15720 

8.8338 

3.3322 

-1.3786 

1.3640 

1.2860 

-0.18934 
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PREE-ELECTRON  MODM..  II 


APPEI^IX 

ILLUSTRATION  OP  THE  METHOD;  NAPHTHACENE 


Matrix  Formulation 

The  naphthaoene  molecule  has  symmetry}  expressed  in  terms  of  cartesian  coor- 
dinates the  symmetry  species  are  SS,  SA,  AS,  AA  (where  the  letters  refer  to  the  behav- 
ior (symmetric  or  antisymmetric)  upon  reflection  at  XY  and  YZ  planes,  respectively). 
See  Pig.  2. 

In  the  case  of  SS  symmetry 


«(1)  n «(4)  - «(5)  - 0(8)  , 

0(2)  - 0(3)  o $(6)  = 0(7)  , 

0(9)  - 0(10)  - 0(11)  a 0(12)  , 

0(1, 12). a 0(4,11)  a 0(5,10)  a 0(8,9)  , 

0(9,12)  a 0(10,11)  , 


(A.l) 


(where '0(r,8 ) designates  the  value  of  the  wavefunotlon  at  the  Joint  atom  located  be- 
tween atoms  r and  s_,  and  Is  to  be  distinguished  from  0(r— t),  the  value  of  the  wave- 
funotlon at  bond  r-t).  Hence,  there  are  only  five  Independent  variables.  If  we  write 
the  secular  equations  [I,  Eq.  (2.44),  (2.45),  (2.46)],  keeping  the  relations  of  (A.l) 
in  mind,  we  arrive  at 


(-3cos(e+l)0(9,12)+  20(12) 

0(9,12)-2oos):0(l2)+  0(1,12) 

0(12)+(1-3oo8k)0(1,12)+  0(1) 

0(1,12)-2oobk0(1)+ 

0(l)+(l-2cos/c)0(2) 


(A. 2) 


By  similar  considerations,  we  can  write  for  AS  symmetry 


-(1+3oo8jc)0(9,12)+  20(12) 

0(9,12)-2oosk0(12)+  0(1,12) 

0(12)+(1-3oosk)0(1,12)+  0(1) 

0(l,12)-2cosic0(l)+  0(2) 

0(l)-(l+2oos/e)0(2) 


(A. 3) 


For  SA  syramettTr, 
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*(9,12)  = *(10,11)  = 0 

-2oo0k*(12)+  *(1.12)  a 0 

*(12)+(1-3c08k)*(1  ,12)+  *(1)  a 0 ■. 

*(1,12)-2co8K*(1)+  *(2)  a 0 

*(1)+(1-2cobk)*(2)  a 0 . 

For  AA  8ytnmatry, 

*(9,12)  a *(10,11)  a 0 

-2oos)e*(12)+  *(1,12)  » o 

*(12)-(1+3oo8/c)«(1,12)+  *(1)  a 0 

*(l,12)-2oos(c*(l)+  (*(2)  a 0 

*(1)-(1+2oo8k:)*(2)  a 0 , 


(A. 4) 


(A. 5) 


The  Eigenvalues 

From  this  point  on  we  will  discuss  only  the  symmetry  SS.  From  (A. 2)  we  obtain  the 
following  secular  equation* 


(l~3oos/e) 

1 


2 

-2cobk 

■1 


(l-3oos/e) 

1 


1 

-2oosk  I 
1 (1-2oo8k) 


(A. 6) 


giving 


-72(oosic)^  + 84(oosk)^  + 34(cobk)^  - 48(cos)e)®  - 3cosk  + 5 a o 


or 


(oobk-1)[72(oo8k)^  - 12(oo8k)^  - 46(oosk)®  + 2oobk  + 5]  a o. 


(In  the  absence  of  free  end  points,  the  solution  cobk  a i oan  be  anticipated  for  S3 
symmetry.)  The  solutions  of  (A. 6)  are  given  in  Table  I for  the  S3  orbitals.  The  neg- 
atives of  these  roots  are  the  roots  of  the  AS  symmetry  secular  equation. 


Wave function  Obtained  vIa  the  Eigenvectors 

The  secular  equation  (A. 2)  may  now  be  solved,  for  each  eigenvalue,  to  give  the 
relative  values  of  the  eigenvectors.  Application  of  the  normal iaat ion  condition  [I, 

Eq,  (2.34)]  fixes  their  absolute  values. 

Thus 
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PREE-ELECTRON  MODEL.  II 
<t>{2)  = (-1+2oosk)"^0(1)  , 

0(1,12)  = [2coaK-(-l+2cosK)“^]0(l)  , 

0(12)  B {(-1+5oo8k) (2oosk-(-1+2cosk)“^]-1}0(1)  , 

0(9»12)  = 2{(-l+3CO3/c)t2COS/C-(-l+2COSK)"^]-l}(-l+3CO3K)"^0(l)  , , 


(A. 7) 


40®(1)  + 40^(2)  + 40^(12)  + 4(3/2)0^(l,12)  + 2(3/2)02(9,12)  » 1 . (A. 8) 

For  example,  for  the  eigenvalue  eoo/c  a i 


0(r)  a 0(1)  for  all  r, 

02(1)  a 1/21  a 0.04762. 

For  the  eigenvalue,  oobk  a 0.79002, 

«(1)  - «(1) 

0(2)  a 1.72400(1)  , 

0(1,12)  a -0.14400(1)  , 

0(12)  . -1.19730(1)  , 

0(9,12)  a -1.74780(1)  , 

and  from  (A. 8) 

02(1)  a -235  , 0(1)  - *0.1799  . 

If  we  arbitrarily  aeleot  0(1)  to  be  positive,  all  the  other  eigenvectors  are  uniquely 
determined.  We  may  now  find  the  value  of  the  wavef unction  at  the  mld-polnt  of  a bond 
from  Eq.  (3a).  We  can  find  a^  and  from  the  considerations  of  I,  Eq.  (2.21)-(2.30) 
and  the  accompanying  text;  of  and  tanfig  (the  quantities  listed  In  Tables  II  and 
III)  from  Eq.  (3b)  and  (3o),  and  a^  and  from  these.  Since  tt/2  63  > -V2,  we 
obtain,  for  example,  for  cobk  a 0.79002, 


6 

6 

6 

6 


1 

2 

3 

4 


0 , 

0.8296  , 

0 , 

0 , 


(A. 9a) 
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In  radians.  Hence,  0.1799  = 0(1)  = cos(3(c/2)  = 0.5488,  etc. . so  that  one  finds 


- 0.3278  , 

Sg  = “O.3189  t 
a^  o -0.0274  , 
aj^  = -0.3323  • i 


(A. 9b) 


Wave function  Obtained  via  the  Joint  Conditions 

This  method  works  without  excessive  labor  If  the  molecule  is  small  and/or  there 
Is  enough  symmetry  to  supplement  It.  Consider  the  case  In  question  (S3  symmetry). 
Here  6^^  ■ 6^  a ■ 0,  since  we  have  chosen  our  branch  origins  to  lie  on  the  symmetry 
axis.  (If  the  axis  had  been  an  A axis,  we  would  have  6^  ■ tt/2.)  There  remains  6g 
to  be  determined.  We  write,  e.£. , at  the  Joint  atom  (9»12)  the  boundary  conditions 


:2tan(-/e+6g)+tan(K/2)  » 0 , (A.lOa) 

ag0os(-/t+ag)  a a^oos()c/2)  . (A.  10b) 


By  use  of 

tan(iex3+ag)  a (tan/cxg+tana2)/(l-tanxxg  tan6g)  , (A. 11) 

Eq.  (A.lOa)  is  transformed  from  a transcendental  equation  in  6g  into  an  algebraic 
equation  In  tandg,  which  la  readily  solved  to  yield  the  value  in  (A. 9a).  By  virtue 
of  (A. 11),  the  system  of  Joint  conditions  Is,  In  all  cases,  equivalent  to  a system  of 
algebraic  eqmitlon. 

The  values  of  the  when  substituted  into  the  equations  of  the  type  (A.lOb), 
readily  determine  the  ratios  of  the  a^'s.  The  normalization  condition  I,  Eq.  (2.7) » 
then  gives  the  absolute  values  listed  In  (A. 9b). 

The  value  of  the  wavefunctlon  at  any  point  Is  now  easily  found,  e_.fi,., 

at  atom  12:  0(12)  » agoosdg  , 

at  atom  9,12:  0(9,12)  «>  agOOB(-)c+6g)  ■ a|^eos(K/2)  , 

at  bond  9,12-12:  0(9,12-12)  - agC03(-)c/2+6g) . 

Use  of  Tables  I,  II,  III  for  the  Construction  of  the  Wave functions 

The  fig's  may  be  obtained  from  the  tanfig's  listed  in  Table  III.  The  signs  of  the 

fig's  are  the  same  as  the  signs  of  the  tanfig's  [See  I,  Eq.  (2.11)].  The  magnitudes  of 

o 

the  Sg's  are  the  square  roots  of  the  a^  '3  listed  in  Table  II.  We  arbitrarily  choose 
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one  of  the  a2(sRy  to  be  positive.  To  find  the  sign  of  Sg,  we  write  the  continuity 
condition  I,  Eq.  (1.32)  or  (2.4),  at  atom  (1,12) « 


&^oob5k/2  = agOos(K+6g)  . 

The  fllgna  of  the  cosines  determine  the  sign  of  the  ag.  We  may  then  oonsiv.er  (A. 10b) 
to  find  the  sign  of  a^^,  etc. . . . 

Convenient  checks  on  the  calculations  are  provided  by  I,  Eqs.  (2.43),  (2.89), 
and  (2.90). 


